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Using dispersion relations, we derive the complete virtual QED contributions to Bhabha scattering 
due to vacuum polarization effects. We apply our result to hadronic corrections and to heavy lepton 
and top quark loop insertions. We give the first complete estimate of their net numerical effects 
for both small and large angle scattering at typical beam energies of meson factories, LEP, and the 
ILC. With a typical amount of 1-3 per mille they are of relevance for precision experiments. 
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The determination of the luminosity at lepton and 
hadron colliders is a necessary task, since in many cases 
the normalization of the measured cross sections is an 
observable of direct phenomenological interest. In prac- 
tice, this task can only be solved by selecting a particu- 
lar reference process, which is expected to generate large 
statistics, be as free as possible of systematic ambigui- 
ties and predicted by the theory to suitable accuracy. As 
far as lepton colliders are concerned, the above criteria 
are fulfilled by Bhabha scattering, i.e. the e+e" e+e" 
process, where a precision under the per mille level can be 
achieved on both the theory and the experimental sides 

In the last few years, there has been major progress 
in the evaluation of the corrections at the next-to-ncxt- 
to-leading order accuracy. In fact, the two- loop QED 
corrections were first evaluated in the massless case in 

. The photonic corrections to massive Bhabha scatter- 
ing with enhancing powers of ln(s/mg) were soon derived 
from that Q . The missing constant terni@] plus the cor- 
rections with electron loop insertions f?, 'sffoUowed later. 
Recently, the heavy fermion (or Nf — 2) corrections were 
derived in the limit ml « m? « s, |t|, |m| [1, @|, where 
m is the mass of the heavy fermion, and soon after also 
for arbitrary m, with ml « m^, s, |m| 

[ia[iii,[ii. 

In this letter, we present the last missing part of the 
virtual corrections, the hadronic ones. 

The three classes of two-loop diagrams that we con- 
sider are shown in Fig. [1] They all may be evaluated by 
dispersion integrals, after replacing the vacuum polariza- 
tion insertion Ilhadl'?^) to the photon propagator |13)], 
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by the once-subtracted dispersion integral 



nhad(<z') = 



TT 



4M? 



dz Imnhad(z) 
z q"^ — z + iS 



(2) 




FIG. 1: Three classes of two-loop virtual hadronic Bhabha 
diagrams, (a) and (b) represent hadronic irreducible and 
reducible vertex diagrams, (c) irreducible hadronic box di- 
agrams. 



ratio i?had, 
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Finally, one relates Im Ilhad to the hadronic cross-section 



measured experimentally in the low-energy region and 
around hadronic resonances, and given by the perturba- 
tivc QCD prediction elsewhere. For heavy fermion inser- 
tions, we have in Eq. fS]) instead of i?had(-z): 

Rf{z) = Q}Cf{l + 2m}/z)^l - Amj/z, (4) 

with electric charge Qf and color factor C/, to leading 
order, which is sufficient for practical purposes. For lep- 
tonic or top quark intermediate states the dispersion rela- 
tion approach is just an efficient technique of evaluation, 
but it becomes essential for light quark loops. In the case 
of Bhabha scattering, this method was first used some 
time ago for one-loop propagator insertions It was 
also applied to two-loop irreducible vertex (plus soft real 
pair) corrections (3]. Here, we derive semi-analytical 
cross section formulae for the so far unknown hadronic 
box diagrams (see Fig.[T}:). 

After replacing the photonic self energy insertion in a 
two-loop diagram by the integral ^ and subsequently 
interchanging loop and dispersion integrations, one ar- 
rives at integrals given by the convolution of the hadronic 
data with kernel functions K{z), where the latter repre- 
sent massive one-loop Feynman diagrams. Details of the 



2 



evaluation will be presented elsewhere. The one-loop self- 
energy kernel is trivial, -ftrsE(-2) = 1/(9^—2)- Whereas the 
two- loop vertex kernel K„ can be found in Eq. (5) of [31 j 
the cross section corrections due to the double boxes of 
Fig . IJb depend on three such kernels Ka{z),a = A,B,C 
[12| [ISj . Notice that, unlike the vertex kernel, the box 
kernels are infrared divergent, but, analogously to the 
one-loop box, they have no singularity in the electron 
mass. The net cross section contribution from the eight 
double box diagrams is still infrared divergent. These 
boxes become, as well as the reducible diagrams with 
one-loop vertices and boxes, infrared finite after adding 
real soft photon emission. The anatomy of that is nicely 
detailed in Section 2.2 of [3]. In order to construct an 
infrared-finite quantity, we combine: (i) Born diagrams 
interfering with the two-loop box diagrams (Fig. [!};) and 
with reducible vertices (Fig.[T]D), (ii) Born diagrams with 
a one-loop vacuum polarization function interfering with 



single one-loop boxes and vertices, and finally (iii) real 
single-photon emission with a one-loop vacuum polariza- 
tion [l9|. The resulting cross-section becomes: 
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with c = a4/(7r2s) and i?h(s) = d{s - 4M^) i?had(s) 
Further, 
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We use the abbreviation c(s,<) = [(s + + (2s^ + 

2si + t^)/i]/3. The real soft photon emission is cut at a 
maximal photon energy E^'^^ = uj. After combining with 
real hard photon emission from a Monte Carlo program, 
the dependence on this parameter disappears. This is 
simulated here, as usually, by setting formally 2ijj/^ — 
1, which switches off the corresponding terms. 

Of course, we can get, from Eqn. ([5]), the Nf = 2 
contributions from heavy fermion loop insertions, with 
the replacements 4A/^ ^m^, and R\iaA{z) Rf{z), 
see Eqn. g]) [2^. 

We will now discuss the numerical net effects arising 
from the Nf = 2 vertex plus box diagrams (i.e. excluding 
the pure running coupling effects): 
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with da/dfl from Eqn. ([5|). The expression for the irre- 
ducible vertex term da^/dQ derives directly from 1^, 16| . 
The da2/dQ is normalized to the pure photonic Bhabha 
Born cross section dao/dQ: 
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One has to compare the da2/dQ with the anticipated 
experimental accuracies varying from few tenths of per 
mille (at small angles e.g. at LEP or ILC/GigaZ) to few 
per mille at large angles (at meson factories, but also at 
LEP and the ILC) [l|, l3il3|]. It also compares to the pure 
photonic two- loop corrections which amount to few 
per mille and are relevant at all energies. The large an- 
gle region at y/s = 1 GeV is studied e.g. at KLOE at 
DA$NE (Fig. [2^). Here it is also important that our 
formulae are valid for an arbitrary ratio s/mj: or s/rn^, 
as long as ml is small. The hadronic corrections are 
very small, the net fermionic ones get about 1 per mille. 
At y/s = 10 GeV (Fig. Wp), the hadronic corrections 
dominate and the net corrections amount to about 3 per 
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FIG. 2: Two-loop Nj = 2 vertex and box corrections da2 
to Bhablia scattering in units of 10~'^do-o at meson factories, 
y/s = 1 GeV (a) and ^/s = 10 GeV (b). 



mille. We also show the small angle region at y/s = 100 
GeV, see Fig. [Sb, which is important for the luminos- 
ity determination at LEP and the GigaZ option of the 
ILC. Largest are the hadronic terms, and the net effect 
amounts to more than 1 per mille. Finally, at yfs = 500 
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FIG. 3: Same as Fig. [1 but ILC energies of ^ = 100 GeV 
(GigaZ option) and ^ = 500 GeV. 



GeV, the net flavor terms amount to about 8 per mille 
and are half of the photonic effect. In [l^, we dis- 
cussed for leptons the individual terms from self-energies, 
irreducible vertices, and the infrared sensitive part with 
boxes separately. For the electron and heavy lepton cor- 
rections, there is also complete agreement with the results 
of a non-dispersive approach . As a further check on 
the consistency of the calculation we note that the cor- 
rections show decoupling when the masses TOj or be- 
come large compared to the kinematics. In fact, the top 
quark and tau lepton contributions decouple nearly 
everywhere; they are included in the total non-photonic 
effects shown in the figures. 

Naturally, the predictions depend on i?had, i-e. on ex- 
perimental data. Although i?had is being used in many 
problems (like the muon anomalous magnetic moment for 
example), no recent fit is publicly available. So, we were 
forced to use an older version of i?had, as it was used in 
(lit , which was available by contacting the author 17 1. It 



is expected that current hadronic data would not induce 
changes larger than about 10% in our analysis. Because 
the size of the hadronic effects here is small enough so 
that only the leading digit truly matters, this is a by far 
sufficient accuracy. 



Summarizing, the net virtual Nf = 2 corrections to 
Bhabha scattering amount typically to 0.1%. They have 
to be taken into account for precision studies, and a pack- 
age for the evaluation has been made public at [3]. For 
the electron and heavy lepton contributions, we agree 
with other computations. The newly evaluated hadronic 
contributions are of comparable size. These virtual cor- 
rections have to be combined yet with cut-dependent real 
fermion pair or hadron emissions. This has to be done 
with MC generators. 

We would like to thank B. Knichl and H. Burkhardt 
for help concerning i?had and A. Arbuzov, H. Czyz, S.-O. 
Moch, and K. Monig for discussions. Work supported 
in part by GGI and INFN in Florence, by SFB/TRR 9 
of DFG, by the Sofja Kovalevskaja Programme of the 
Alexander von Humboldt Foundation, and by MRTN- 
CT-2006-035505 "HEPTOOLS" and MRTN-CT-2006- 
035482 "FLAVIAnet". 



[2: 



[9 
[10 



[11 
[12 

[13 
[14 
[15 

[16 

[17 

[18 
[19 

[20] 



K. Monig, "Bhabha scattering at the ILC", 
Bhabha Workshop, Karlsruhe, April 2005, 



talk at 



http: / /sfb-tr9. particle.uni-karlsruhe.de / 1 



S. Jadach, "Theoretical calculations for LEP luminosity 
measurements", talk at Bhabha Workshop, Karlsruhe, 
April 2005, http://sfb-tr9.particle.uni-karlsruhe.de/ 
Balossini et al, Acta Phys. Polon. B38 (2007) 3441. 
Z. Bern, L. Dixon, and A. Ghinculov, Phys. Rev. D63 
(2001) 053007. 

N. Glover, B. Tausk, and J. van der Bij, Phys. Lett. B516 
(2001) 33. 

A. Penin, Phys. Rev. Lett. 95 (2005) 010408. 

R. Bonciani, A. Ferroglia, P. Mastrolia, E. Remiddi, and 
J. van der Bij, Nud. Phys. B716 (2005) 280. 
S. Actis, M. Czakon, J. Gluza, and T. Riemann, Nucl. 
Phys. B786 (2007) 26. 

T. Becher and K. Melnikov, JHEP 06 (2007) 084. 
S. Actis, M. Czakon, J. Gluza, and T. Riemann, XXXI 
Conference "Matter to the Deepest", Ustron, Poland, 5- 
11 Sep 2007, |http: / /prac.us.edu.pi7~us2'007/talks.htni| 
R. Bonciani, A. Ferroglia, and A. Penin, 
arXiv:0710.4775V3 [hep-ph]. 

S. Actis, M. Czakon, J. Gluza, and T. Riemann, Acta 
Phys. Polon. B38 (2007) 3517. 

N. Cabibbo and R. Gatto, Phys. Rev. 124 (1961) 1577. 
F. Berends and G. Komen, Phys. Lett. B63 (1976) 432. 

B. Kniehl, M. Krawczyk, J. Kiihn, and R. Stuart, Phys. 
Lett. B209 (1988) 337. 

DESY, webpage 

http: / /www-zeuthen. desy.de / theory / research /bhabha/ 1 
H. Burkhardt, TASSO-NOTE-192 (1981), and Fortran 
program repi.f (1986). 

Kernels and other functions may be downloaded at [l6l |. 
The infrared safe Nf = 2 irreducible vertices (see Fig.[T^) 
and pure self energy diagrams are not included here. 
In this case, the kernel integrations in Eqn. ([SJ may be 
easily performed explicitely, leading to analytical expres- 
sions with Euler trilogarithms or simpler functions. 



